Dynamical Creation of Fractionalized Vortices and Vortex Lattices 
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We investigate dynamic creation of fractionalized half-quantum vortices in Bose-Einstein conden- 
sates of sodium atoms. Our simulations show that both individual half-quantum vortices and vortex 
lattices can be created in rotating optical traps when additional pulsed magnetic trapping potentials 
are applied. We also find that a distinct periodically modulated spin-density-wave spatial structure 
is always embedded in square half-quantum vortex lattices. This structure can be conveniently 
probed by taking absorption images of ballistically expanding cold atoms in a Stern-Gerlach field. 
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Topological excitations such as quantized vortices have 
been fascinating for quite a few decades and recently have 
also been thoroughly studied in Bose-Einstein conden- 
sates (BECs) of ultracold atoms 0, H i, i, H, @, 0, & % ■ 
For vortices in single component BECs, the circulation of 
supercurrent velocity v s along a closed curve T around a 
vortex line (defined as C — J r dl- v s ) is quantized in units 
of 2irh/m (to is the atomic mass), with C = ±1,±2, ... 
as a consequence of analyticity of single- valued wavefunc- 
tions of coherent quantum states. Furthermore, only a 
vortex with circulation C — ±1 or an elementary vortex is 
energetically stable. A secondary vortex with circulation 
C = ±2, ±3, ... spontaneously splits into a few elementary 
ones which interact via long range repulsive potentials. 

A configuration with its circulation smaller than the 
elementary value (C = 1) has to be described by a sin- 
gular wavefunction and it always turns out to be ener- 
getically catastrophic. A most obvious example is a two- 
dimensional configuration, where the condensate phase 
angle <&(r, 9) rotates slowly and uniformly by 180° in the 
r — plane around a vortex center but jumps from tt to 
27T when the polar angle is equal to 27r. The 7r-phase 
jump here effectively induces a singular cut in the wave- 
function. The corresponding circulating velocity field is 
simply v s (r, 6) — H/(2mr)eg, leading to C = 1/2 that is 
one-half of an elementary value. The energy of a cut per 
unit length along 9 = 2ir line where the phase jumps is 
finite and therefore the overall energy of a cut in an in- 
dividual fractionalized vortex scales as L, L is the size of 
system, while the energy for an integer vortex only scales 
as a logarithmic function of L. Consequently, a cut that 
connects two singular fractionalized vortices mediates a 
linear long rang attractive potential that confines all frac- 
tionalized excitations. So in a single component conden- 
sate, vortices with C — ±1 are fundamental ones which 
do not further split into smaller constituent elements as 
a result of confinement of fractionalized vortices. 

Hyperfine-spin degrees of freedom can drastically 
change the above arguments about elementary vortices. 
In condensates of sodium ( 23 Na) or rubidium ( 87 Rb) 



atoms in optical traps [Io| , hyperfine spins of cold atoms 
are correlated because of condensation. A pure spin de- 
fect (or a spin disclination in the case of 23 Na) where 
spins of cold atoms slowly rotate but no supercurrents 
flow, can carry a cut, i.e. a line along which a 7r-phase 
jump occurs as a result of Berry's phases induced by spin 
rotations Such a spin defect can then terminate a 



cut emitted from a singular half-quantum vortex (HQV) 
configuration, which consequently leads to a linear con- 
fining potential between the spin defect and HQV. For 
instance, a HQV with C = 1/2 confined to a spin defect 
does exist as a fundamental excitation in spin nematic 
condensates 12, llij . In this Letter we will present the 



simulation of dynamical creation of fractionalized HQ Vs 
in a rotating BEC of sodium atoms, and formulate an 
experimental procedure for the realization of such exotic 
topological excitations. 

Magnetic properties and energetics of half vortices Un- 
like conventional integer vortices, HQVs have very rich 
magnetic structures. We will demonstrate that far away 
from its core a HQ V have vanishing local spin densities 
but is accompanied by slowly rotating spin quadrupole 
moments; within the core, spin densities are nonzero. 
The Hamiltonian for interacting sodium atoms is 
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Here t^ai^a), a — x,y, z are creation (annihilation) op- 
erators for sodium atoms in hyperfine states \a); they 
arc defined as linear superpositions of creation opera- 
tors for three spin-one states, |1,to^), mp — 0, ±1. 
V4 = (Vl-^i)/V2, vb\ = (v4+^! 1 )/tV2andV! = 4- 
Co, 2 are interaction parameters that depend on two-body 
s-wave scattering lengths ao,2 for total spin 0,2: Co = 
4irh 2 (ao + 2a2)/3m and C2 = 47r/i 2 (a2 — &o)/3to. For con- 
densates of 23 Na atoms, ao — 50 as and — 55 as (as 
is the Bohr radius). S a = — ze Q/ g 7 5'^5' 7 and p = 4' a ' t Pa 
are local spin-density and density operators. 

Generally, a condensate wavefunction ^> a {= (^l)) for 
spin-one atoms is a complex vector. For sodium atoms, 
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interactions favor states with a zero spin density; this 
leads to a spin nematic ground state which does not break 
the time reversal symmetry and VP = exjp(i&)^/pn; here 
n is a unit vector with three components n a , and p is 
the number density of sodium atoms. This spin nematic 
state is invariant under an inversion of n and a 7r-phase 
shift (i.e. n — * — n and d> — * $ + 7r). Although local 
spin densities (S(r)) in a nematic state vanish, a nematic 
condensate carries a spin quadrupole moment defined as 
Qap{ Y ) — (SqS^) — i<5 Q/ 3(S 2 ), which can be calculated 
and is specified by the nematic unit director n introduced 
above: Q a p = p{p. a Hj3 - \8 a p). 

Around a HQV centered at the origin and oriented 
along the z-direction, nematic directors lying in a per- 
pendicular r — 9 plane rotate by 180° forming a spin- 
disclination; the corresponding condensate wavefunc- 
tion far way from the vortex core is ^(0, r — > oo) = 
exp(i9/2)y/pn(9), n x =cos(9/2), n y =sm(9/2) and n z =0. 
The 180° rotation of nematic director n(8) around the 
vortex illustrates that a 7r-spin disclination where spin 
quadrupole moments Q a/ 3 (9) slowly rotate is indeed con- 
fined to this HQV. In the Zeeman basis of |1, mp) states, 
772 p = 0, ±1, the above vortex state is equivalent to ip\ = 
exp(i6»)/(r)^p/V2, ip-i = -g{r)^/~p/\/2 and Vo = 0; far- 
away from the vortex core, f(r — > oo) = g(r — > oo) = 1. 
The core structure can be studied by numerically solving 
the multiple-component Gross-Pitaevskii (GP) equation 

h 2^;V 2 + (c + c 2 )p±i + (co - c 2 )p T i}ifj±i = 0, (2) 

here p mF = \ip mF \ . The corresponding boundary condi- 
tions when r — > oo are set by the asymptotic behaviors of 
a HQV fax away from the core as discussed before Eq.([2]). 
Notice that the last term in Eq. ^ indicates that mutual 
interactions between |1,±1) atoms induced by scattering 
are repulsive since cq — c 2 is positive. 

Our results show that within the core, |1,1) atoms 
state are completely depleted while |1, — 1) atoms are 
not. This is because supercurrents are only present in 
|1, 1) component; in fact the density of |1, —1) atoms has 
an additional small bump at the center of core, to fur- 
ther take advantage of the depletion of |1,1) atoms in 
the same region to minimize the overall repulsive inter- 
actions between |1, ±1) atoms. So the core has a nonzero 
local spin density (S z (r))(= |V>i| 2 — l^-il 2 ) with excess 
atoms at |1, — 1) state and we define it as a minus HQV. 
Similarly, one can construct a plus HQV with an identical 
vorticity (i.e. V x v s ) but with excess |1, 1) atoms in its 
core (See FigQJa)). Generally, HQVs have distinct spa- 
tial magnetic structures: a HQV core carries excess spins 
while far away from the core spin quadrupole moments 
slowly rotate around the vortex. 

Interactions between two HQVs very much depend on 
species of HQVs involved. When both are plus ones (or 
minus), the corresponding interaction potential V P p( mm ) 
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FIG. 1: (color online), (a) Density (p±i = IV'iil' 2 ) an d spin- 
density (pi — p-i) profiles of an individual plus half-quantum 
vortex (HQV) centered at r = 0. (b) Interaction potentials (in 
units of chemical potential p) between two HQ Vs as a function 
of separation distance d. V vm is the potential between a plus 
HQ V and a minus HQ V; V pp is the strong repulsive potential 
between two plus HQVs. Inset is for spin densities in a pair of 
plus-minus HQVs separated at d=30£, £ is the healing length. 

is a logarithmic long range one due to interference be- 
tween coherent supercurrents. However, the interaction 
between a plus and a minus HQ V V pm is repulsive and 
short ranged only extending over a scale of vortex cores. 
In this case, supercurrents flow in different components 
and they don't interfere; the short range interaction is en- 
tirely due to inter-component interactions between |1, 1) 
and |1,— 1) atoms. Indeed, the amplitude of potential 
Vp m is proportional to, when it is small, cq — c 2 which 
characterizes mutual interactions between |1, ±1) atoms. 
An integer vortex or a pair of plus-minus HQ Vs centered 
at a same point therefore is unstable and further frac- 
tionalizes into elementary HQVs. In Fig[T] we summarize 
results of an individual HQV and two HQVs. 

Dynamical creation of HQ Vs in rotating BECs To dy- 
namically create HQVs, we numerically solve the time- 
dependent coupled GP equations of spin-1 BEC 

(* ~ l)^^- = h^V 2 + V tr - p T A - QL Z + c oP 
at 2m 

+ c 2 (p±! + po - p T i) + W±]ip± 1 +c 2 ip 2 'iijz F i, 

( l ~ ^ h ^T = h|~ v2 + V tr -P- VL Z + c p 
at 2m 

+ c 2 (px +p-i)]ip +202^1^-1^0, (3) 

where p = J2 m P™f i s the total condensate density, 
Vtr (r) is a spin- independent confining potential of an op- 
tical trap, and W±(r) are pulsed magnetic trapping po- 
tentials which we further apply in order to create HQVs. 
p and A are the Lagrange multipliers used to preserve the 
total number and magnetization of atoms respectively; 7 
is a phenomenological damping parameter which is nec- 
essary for studies of quasi-stationary states 0]. 

We restrict ourselves to a cigar-shaped potential with 
the aspect ratio A = ujj_/uj z ~ 14 which was also used 



3 



0.02 i 




ViDV, 



-5 





5 














-5 





5 



FIG. 2: Creation of a half-quantum vortex. Density pro- 
files of |V+i| 2 , |^-i| 2 , \^+i\ 2 +\^-i\ 2 , and spin-density profile 
| 2 -|?/)_i | 2 (bold line) are shown. The rotating frequency 
is suddenly decreased from an initial value = 0.65tJ_L to 
Q= 0.3u± at t — S00 ms. The bottom panel shows that a 
single half-quantum vortex is formed at i=1600 ms after the 
magnetic trapping potential has been adiabatically switched 
off. 



in early experiments Q ■ We consider a two-dimensional 
cylindrical trap which is characterized by two dimension- 
less parameters C = ** N l£+ 2a *) and C 2 = ^ N ^~^ , 
with L z the size of system along the z-axis and N — 
3 x 10 6 the total number of sodium atoms. When com- 
bined with a nonaxisymmetric dipole potential that can 
be created using stirring laser beams 0, 0] , the optical 
trapping potential in a rotating frame is given by Vt r (r) = 
mw}{(l + e)x 2 + (1 - e)y 2 }/2. Here w± = 2tt x 250 Hz, 
and anisotropic parameter is set to be e = 0.025. We 
also include an additional magnetic trapping potential: 
W±(r) = : f(3muj 2 _(x 2 + y 2 )/2, which could be realized in 
an Ioffe-Pritchard trap via a Zeeman splitting mpgFtiBB 
with the Lande factor gp=—l for sodium atoms. 

We start our numerical simulations with an initial state 
where |1,±1) are equally populated. Experimentally, it 
was demonstrated that when the bias field is small and 
the gradient field along the z-axis of the trap is almost 
canceled, such a state can be prepared and the |1,±1) 
components are completely miscible (however, with the 
immiscibility between |1, ±1) and |1, 0) components) fioj ]. 
We then study the time evolution of this initial state 
using the Crank-Nicolson implicit scheme pj ■ The unit of 
length is ah—y/ h/2muj±— 0.48 /im and the period of the 
trap is uj 1 = 4 ms; the interaction parameters are Co=500 
and C2=450 and the damping rate is 7 = 0.03. Further, 
we include symmetry breaking effects in our simulations 
by allowing the trap center to randomly jump within a 
region [—(5, <5]x[— S, S] (S — O.OOlh, where h is the grid size), 
which is crucial for vortices to enter the condensate one 
by one rather than in opposing pairs @, S] • 

However, without additional pulsed magnetic poten- 
tials, one can show that dynamic instabilities for creation 
of integer vortices in rotating BECs occur almost at same 
frequencies as for HQVs and a triangular integer- vortex 
lattice is formed (see FigJ3Ja)). This integer- vortex lat- 



tice is locally stable with respect to the non-magnetic 
perturbations, by applying an additional optical trap- 
ping potential with an oscillating trapping frequency to 
effectively shake them, indicating their metastability. For 
this reason, a time-dependent magnetic trapping poten- 
tial with harmonic form is applied; and we find that 
when a pulsed magnetic field with (3 > 0.005 is applied, 
HQV lattices could be formed. Here we set (3 = 0.1, 
which is suitable for both generating a single HQ V and 
demonstrating the dynamical evolution of HQ V lattices 
formation. After the magnetic trapping potentials are 
on, 1 1,1) component further spreads to the edge, while 
1 1,-1) component remains at the center of the trap and 
surfaces of equally populated |1, ±1) components become 
mismatched with two different Thomas-Fermi radii. 

First, we dynamically create a single half-quantum vor- 
tex in condensates, which can be used for the study of 
dynamics of a HQV. We switch on abruptly a rotating 
drive with fi = 0.65wj_ and the trap anisotropy e is in- 
creased rapidly from zero to its final value 0.025 in 20 
ms. At t — 800 ms, only one vortex in |1, — 1) component 
appears. Afterwards, We decrease to 0.3lu± suddenly, 
and switch off the magnetic trapping potential adiabat- 
ically within 200 ms. We then find a stable single HQV 
formed at t = 1600 ms as shown in Fig|2] The 0.3ui± fre- 
quency used after t = 1600 ms is within the stable region 



estimated earlier 14j and our simulations of dynamics 
are consistent with the energetic analysis. 

Fractionalized-vortex lattices can be created in a sim- 
ilar setup. The main experimental procedure and re- 
sults of our simulations for creation of HQV lattices 
are presented in Figl3] After a rotation with frequency 
f2 = 0.7cjj_ starts abruptly and anisotropy e is set to 
its final value 0.025, we can see that the cloud is initially 
elongated and at the same time rotates with the trap. At 
about 150 ms, surface ripples due to quadrupole excita- 
tions occur in |1, 1) component of the condensate, while 
no surface oscillations appear on the surface of |1, — 1) 
component. At t = 240 ms, we find that the density 
profile of |1, 1) component is along a short-axis while the 
1 1,-1) component is along a long-axis due to repulsive 
interactions between two components, and the surface 
of 1 1,-1) component is not always buried in the inner 
region of |1,1). The surfaces of two components oscil- 
late independently and are decoupled dynamically. At 
t = 430 ms, we find that two plus HQVs with excess 
|1,1) atoms inside cores have nucleated at the center. 
Correspondingly, we observe two small regions near the 
center where |1,— 1) atoms are completely depleted and 
the density of |1, 1) atoms remains high. At t = 800 ms, 
two components are phase separated, but the structure of 
HQ Scores remains almost unchanged. Finally, we switch 
off the additional magnetic potential adiabatically within 
200 ms and a HQ V lattice with interlaced square config- 
uration becomes visible. In this structure, to minimize 
strong repulsive interactions V^ p ( mm ) between plus or mi- 
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FIG. 3: (a) Creation of a triangular integer-vortex lattice in 
a rotating optical trap at t — 1600 ms. (b) Creation of half- 
quantum vortex lattices when an additional pulsed magnetic 
trapping potential is applied. Here time evolution of various 
condensate densities is shown. The optical trap rotates at fl — 
0.7lo± with a magnetic trapping potential on until t — 800 
ms; afterwards, the magnetic trap is adiabatically switched 
off within 200 ms. The bottom panel shows the half-quantum 
vortex lattice formation at t = 1600 ms; a square lattice in 
the spin-density profile is clearly visible. 



nus HQ Vs, the vorticity is evenly distributed among plus 
and minus HQVs, or between |1,±1) components. This 
also indicates that spatially each plus HQ V prefers to be 
adjacent to minus HQVs and vice versa to avoid stronger 
interactions Vpp,V mm and to take advantage of relative 
weaker interactions V pm (See FigQ]). To further minimize 
repulsive interactions V pm between nearest neighboring 
vortices, a plus HQV is displaced away from adjacent 
minus HQVs by a maximal distance. Generally because 
of the asymmetry between V pp and V pm , a bipartite vor- 
tex lattice should be favored over frustrated geometries 
such as triangular lattices where a plus HQV could be 
adjacent to another plus H Q V resulting in much stronger 
repulsion. In our simulations of 23 Na atoms in rotat- 
ing traps with cq ~ 30c2, square vortex lattices have al- 
ways been found. Equilibrium energetics of rectangular 



or square lattices were also considered in the quantum 
Hall regime [HI, [l6j], in two-component BECs coupled 
by an external driving field where vortex molecules are 
formed [TtJ , and also observed in condensates of pseudo- 
spin- 1/2 rubidium atoms (l8j . Here we have mainly fo- 
cused on dynamical creation of HQ V lattices confined to 
a spin-density-wave structure at relatively low frequen- 
cies; this structure can be conveniently probed by taking 
absorption images of ballistically expanding cold atoms 
in a Stern-Gerlach field [T3| • 

In conclusion, we have demonstrated a practical setup 
to create fractionalized vortices and vortex lattices in 
BECs of sodium atoms. We found that a square half- 
quantum vortex lattice has a distinct periodically modu- 
lated spin-density-wave spatial structure, due to short 
range repulsive interactions between neighboring half- 
quantum vortices. Our results are of particular signifi- 
cance for creating these excitations in experiments and 
for exploring novel phenomena associated with them. 
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